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In the present work, we study the quantum cosmology description of closed Friedmann- 
Robertson- Walker models in the presence of a positive cosmological constant and a generic 
perfect fluid. We work in the Schutz's variational formalism. If one uses the scale factor and 
its canonically conjugated momentum as the phase space variables that describe the geomet- 
rical sector of these models, one obtains Wheeler-DeWitt equations with operator ordering 
ambiguities. In order to avoid those ambiguities and simplify the quantum treatment of the 
models, we introduce new phase space variables. We explicitly demonstrate that the trans- 
formation leading from the old set of variables to the new one is canonical. In order to show 
that the above canonical transformations simplify the quantum treatment of those models, 
we consider a particular model where the perfect fluid is dust. We solve the Wheeler-DeWitt 
equation numerically using the Crank-Nicholson scheme and determine the time evolution 
of the initial wave function. Finally, we compare the results for the present model with the 
ones for another model where the only difference is the presence of a radiative perfect fluid, 
instead of dust. 
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I. INTRODUCTION 

Since the pioneering work in quantum cosmology due to DeWitt [1], many physicists have worked 
in this theory. The main motivation behind quantum cosmology is a consistent explanation for 
the origin of our Universe. So far, the most appealing explanation is the spontaneous creation 
from nothing In that picture for the origin of the Universe, the Universe is a 

quantum mechanical system with zero size. There is a potential barrier that the Universe may 
tunnel with a well-defined, non-zero probability. If the Universe actually tunnels, it emerges to 
the right of the barrier with a definite size. The application of the creation from nothing idea in 
minisuperspace models has led to several important results. The wave-function of the Universe 
satisfies the Wheeler-DeWitt equation ^j]. Therefore, one needs to specify boundary conditions 
in order to solve that equation and find a unique and well-defined wave function. The motivation 
to obtain a wave-function that represents the creation from nothing has led to the introduction of 
at least three proposals for the boundary conditions for the wave- function of the Universe [7J. The 
inflationary period of the Universe appears very naturally from the creation from nothing idea. 
That is the case because most of the minisuperspace models considered so far have a potential 
that decreases, without a limit, to the right of the barrier. It gives rise to a period of unbounded 
expansion which is interpreted as the inflationary period of the Universe Q]. Also, it was shown by 
several authors that an open inflationary universe may be created from nothing, in theories of a 



single scalar field for generic potentials 



9 



10, 



111 ]. Another important issue is the particle content 
in the Universe originated during the creation from nothing process H . . 

In the present work, we study the quantum cosmology description of Friedmann-Robertson- 
Walker (FRW) models in the presence of a positive cosmological constant (A) and a generic perfect 
fluid. The spatial sections of the models have positive curvature. The perfect fluid has an equation 
of state of the form p = aw, where p and w are, respectively, the fluid pressure and density, 
and a is a constant in the range — 1 < a < 1. In the present quantum cosmology description, 
the perfect fluid is treated by means of the variational formalism developed by Schutz []jj]. The 
treatment of a simpler model which differs from ours by the lack of a cosmological constant was 

n 

first done in [15|. The Schutz's variational formalism, applied to the present model, leads to a 
superhamiltonian which, after the canonical quantization, produces a Wheeler-DeWitt equation 
with operator ordering ambiguities. More precisely, the canonically conjugated momentum to the 
scale factor (p a ) appears in the superhamiltonian coupled to a function of the scale factor (a). 
In the present work we propose a general canonical transformation in order to eliminate these 
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operator ordering ambiguities. As we shall see, the resulting Wheeler-DeWitt equation in the 
new variables will be much simpler than the initial one. In fact, there is already in the literature a 
particular example of such canonical transformation [if]]. There, the authors considered a flat FRW 
geometry coupled to a dust perfect fluid (a = 0) and a cosmological constant. After the application 
of Schutz's variational formalism and the particular canonical transformation, the Wheeler-DeWitt 
equation in the new variables was reduced to a Schrodinger equation for a harmonic oscillator when 
A < and for an inverted harmonic oscillator when A > 0. Although the authors of [lf| state 
that the transformation they used is canonical, they did not prove it. Here, we prove that the 
transformation leading from the scale factor and its canonically conjugated momentum to a new 
set of canonically conjugated variables, such that in the new variables the resulting Wheeler-DeWitt 
equation is free from operator ordering ambiguities, is indeed canonical. In order to prove that the 
transformation is canonical we shall demonstrate that the Dirac and Poisson brackets of the new 
canonical variables are the same. It will be done following the symplectic framework \v\ . Then, we 



16| , on how these transformations simplify 
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shall give a further example, than the one given in Ref. 
the Wheeler-DeWitt equation of these models. We shall consider a closed FRW model coupled to 
a dust perfect fluid and a positive cosmological constant. Finally, we compare the results for the 
present model with the ones for another model already treated in the literature 
difference is the presence of a radiative perfect fluid, instead of dust. 

In the next Section, we write the superhamiltonian constraint in terms of the scale factor, the 
variable associated to the perfect fluid and their conjugated momenta. Then, the superhamiltonian 
term which will lead to operator ordering ambiguities, after the quantization, will be easily seen. 
After that, we introduce the transformations leading to new pairs of canonical variables such 
that the new quantum theory is free from operator ordering ambiguities. Next, we demonstrate 
that the Dirac brackets of the new canonical variables are equal to their corresponding Poisson 
brackets, which confirms that the transformations, leading from the old variables to the new ones, 
are canonical. This will be done following the symplectic framework [171 ]. In Section [Till we give 
an explicit example on how the canonical transformation simplifies the treatment of the model 
at the quantum level. We consider a specific example where the perfect fluid is dust (a = 0). 
We canonically quantize the model, written in terms of the new pair of canonically conjugated 
variables, and obtain the corresponding Wheeler-DeWitt equation. We solve it, numerically. For 
particular values of the dust energy and the cosmological constant, we plot the square modulus of 
the wave- function of the universe as a function of the scale factor. The tunneling process can be 
readily seen from the results. Then, we evaluate the tunneling probability (TP) and its dependence 
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on the dust energy. We obtain that the TP increases with the dust energy for a fixed cosmological 
constant. Therefore, it is more probable that the classical evolution should start with the greatest 
possible value for the dust energy. Finally, we compare the tunneling probability obtained here 
with the one obtained for a model where the perfect fluid is radiative (a = 1/3) [18j]. We find that 
the tunneling probability of the universe with dust is greater than the universe with radiation. 
Finally, in Section [IV] we summarize the main points and results of our paper. 



II. THE CLASSICAL MODEL AND THE CANONICAL TRANSFORMATION 

The closed Friedmann-Robertson- Walker cosmological models are characterized by the scale 
factor a{t) and have the following line element, 



ds 2 = _ N 2 {t)dt 2 + a 2 {t) + r 2^ 2 j (1) 

where dVL 2 is the line element of the two-dimensional sphere with unitary radius, N(t) is the lapse 
function and we are using the natural unit system, where H = c = G = 1. The matter content of 
the model is represented by a perfect fluid with four-velocity U 11 = 5$ in the comoving coordinate 
system used and a positive cosmological constant (A). The total energy- momentum tensor is given 
by 



T^v = (w + p)U^U u -pg^v - Afl^j, , (2) 

where w and p are, respectively, the fluid density and pressure. Here, we assume that p = aw, 
where —1 < a < 1, which is the equation of state for a perfect fluid. 

Einstein's equations for the metric (pQ) and the energy momentum tensor ([2]) are equivalent to 
the Hamilton equations generated by a total Hamiltonian. In order to get this Hamiltonian, the 
symplectic method will be applied in a seminal Lagrangian that describes the dynamics of the 
system, 
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which is an extended formulation of Schutz's Lagrangian, since the a parameter was introduced 
into the model. Calculating the canonical momenta, namely, 
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the Lagrangian density, eq.©, can be rewritten in a first-order form, given by 



C = P a a + P e e + P B B + 



NP n 



12a( 2 " 3Q ) 

Introducing the following transformation 
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the Lagrangian density, eq.©, becomes 



p £ 4/ V 
31/3 ' 

P^e B 



2.3V3 



T, 



£(0) = Pa a + P T r + ^3 a) + 3Na 3a - ANa^ 2+3 ^ - N^, 



= P a a + P T T-V, 
where V = NO, is the symplectic potential, with 



n 



3a 3a + Aa( 2+3a ) + ^. 



12a( 2 " 3a ) a 
Now, we are ready to apply the symplectic method. The symplectic variables 
£^ = (a, P a , T, Pj>, N) and the one- form momenta (A{) are given by 



A a 


= Pa, 


A Pa 


= 0, 


A T 


= Pt, 


A Pt 


= 0, 


A N 


= 0. 
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The corresponding symplectic matrix, denned as 



f, 



dAj dAi 



(12) 



reads 



'0-10 
1 
0-10 
10 

This matrix is singular, so it has a zero-mode, 



(13) 



V 



^=(^0000 1 J . (14) 

Following the symplectic method, if one contracts the zero-mode with the symplectic potential 
gradient one should find zero. Otherwise, a constraint is obtained. Contracting v eq. (|14p with the 
symplectic potential gradient, we get 



u~ = n, (15) 

OXi 

which is a constraint. This constraint will be introduced into the first-order Lagrangian, eq.®, 
through a Lagrange multiplier, /3, 



£« = P a a + P T t + SI/3-V. (16) 

Now, the symplectic variables are Q = (a, P a ,T, Rp, N, (3) and the respective symplectic matrix 
is given by 
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This matrix is singular and, therefore, has a zero-mode, 



an an 



an 



1 1 



(18) 



Contracting this zero-mode with the symplectic potential gradient, we get 



Ki~ — r^. 

OXi 



(19) 



This constraint was obtained before and, following the symplectic method, the system has a gauge 
symmetry. In order to quantize the system, the symmetry must be fixed. This is done in the 
symplectic process introducing the gauge fixing term into the original first-order Lagrangian, eq.([9]), 
through a Lagrange multiplier, r/, 



6 2) = P a d + P T f + St) - V, 



(20) 



where the gauge fixing term is 



Y, = N-a. 



(21) 



(2) 

The symplectic variables are, now, Q = (a, P a ,T, Pj>, N,rf) and the symplectic matrix is 



( -1 

1 

0-10 

1 

1 

y i o o o -l o. 



1 



(22) 



/ 



This is a nonsingular matrix and, according to the symplectic process, its inverse allows us to get 
the Dirac brackets; the only non- vanishing brackets are 



{a,P a } = {T,P T } = l. 
The symplectic potential is identified as being the Hamiltonian, so 



(23) 
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H = - 12 Ji! 3a) " 3a 1+3a + Aa (3+ ^ + P r . (24) 

where P a and are the momenta canonically conjugated to a and T, the latter being the canonical 
variable associated to the fluid [151 ] . This is the total Hamiltonian that generates Hamilton's 
equations equivalent to Einstein's equations for the metric ([TJ and the energy-momentum tensor 
([2]). Since a quantization process will be applied, the first term in the total Hamiltonian (|24p poses 
an operator ordering ambiguity. In order to solve this problem, the following transformation will 
be applied, 



3fl - a)x\ 2/3{1 - a) 
0= ( 1 } ) • (25) 
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Due to this, the Lagrangian density, eq.©, becomes 



3 /3(l-a).\ 2/3W . 2 Q . r /3(l-a)z\ 2a /( 1 - a > , /3(1 - 2( ' ' 
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4 

„-3B 



(26) 



Using the transformation given in eq.© and with the aid of the momenta P x , P e , Pb computed 
from £ (°) (|26p . the first-order Lagrangian will become 

£'(!) = p x ± + p T T - V, (27) 
where the symplectic potential is V = Nil, with 



o _ P* 2 / 3(l-a) ^ 2Q/(1 - a) 

12 ^ 3(l-q)x N 2 /3(l-a) ^ ^ 2 

After that, the symplectic method could be applied. The symplectic variables are now 
'(1) 

£2 = (x,P x ,T, Pr, N) with the corresponding one-form momenta, 
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A — P 

A Px = 0, 

A T = P T , (29) 
A Pt = 0, 
A N = 0. 

(30) 

Computing the corresponding symplectic matrix, using eq. (|12|) . one finds that it has the following 
zero-mode, 



^=(^0000 1 J • (31) 
Contracting this zero-mode with the symplectic potential gradient, a constraint is obtained, namely, 

dV 

*5S " a (32) 

This constraint will be introduced into the first-order Lagrangian (|27p through a Lagrange multi- 
plier, f3, yielding 



£'( 2 ) = p a a + P T f + Slj3 - V. (33) 

'(2) 

The symplectic variables are, now, ^ = (a, P a ,T, Pp, N, p) and the respective symplectic matrix 
is singular with the following zero-mode, 



an an an n i i i (lA) 

WP„ fin ftPrr, / ' V ' 



dP a da dP T 

Contracting this zero-mode with the gradient of the symplectic potential, we get 

dV 

Kl — = n. (35) 

This constraint was obtained before and, in agreement with the symplectic method, the system has 
a gauge symmetry. This symmetry must be fixed, to this end a gauge fixing term (E) is introduced 
into the original first-order Lagrangian (|27p . 
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C = P x x + P T T + Sr) - V, 



(36) 



with 



2 

Now, the symplectic variables are £j = (x, P X ,T, Pp, N, rf) and the corresponding symplectic matrix 
is nonsingular. From the inverse of the symplectic matrix the non-vanishing Dirac brackets 

{x,P x } = {T,P T } = l (38) 

are obtained; the remaining brackets are all null. The symplectic potential is identified as being 
the Hamiltonian, so 



H = Nn = - Veff. + Ft, (39) 

with 



/3C1 - a )x\ ( 2 + 6Q )/ 3 ( 1 - a ) /3fl _ Q W\ 2(l+a)/(l-a) 

Note that the Hamiltonian above is a Schrodinger-like equation with no operator ordering ambi- 
guities. Further, the Dirac brackets are equal to the Poisson brackets, allowing us to conclude that 
the variable transformation, eq. (|25p . is a canonical transformation. 



III. THE QUANTIZATION OF A PARTICULAR MODEL 



In order to show the usefulness of transformation (|25p . let us consider a particular example. 
Consider a model with closed FRW geometry coupled to a dust perfect fluid (a = 0) and a positive 
cosmological constant (A). The total Hamiltonian, eq. (139p . reduces to 



12 V 2 J V 2 J 

The model will be quantized following the Dirac formalism for constrained systems [ijj]. First 
we introduce a wave-function which depends on the canonical variables x and T, 
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* = *(£,T) . (42) 

Then, we impose the appropriate commutators between the operators x and T and their respective 
conjugate momenta P x and Py. Working in the Schrodinger picture, the operators x and T are 
simply multiplication operators, while their conjugate momenta are represented by the differential 
operators 

ft— Pt (43) 

Finally, we demand that the operator corresponding to H annihilate the wave- function \&, which 
leads to Wheeler-DeWitt equation 

(l2a?- 3 (Tj + A (t) )»(».r) = - S :»(x,r), (44) 
where the new variable r = — T has been introduced. 

n 

The operator TL is self-adjoint [20] with respect to the internal product, 



(*,$)=/ #(s,t)*$(x,t), (45) 

«/ 

if the wave functions are restricted to the set of those satisfying either \l/(0, r) = or ^'(O, r) = 0, 
where the prime / means the partial derivative with respect to x. Here, we consider wave functions 
satisfying the former type of boundary condition and we also demand that they vanish when 
x — > oo. 

The Wheeler-DeWitt equation (|44p is a Schrodinger equation for a potential with a barrier. 
We solve it numerically using a finite-difference procedure based on the Crank-Nicholson method 



211 ] . [221 ] and implemented in the program GNU-OCTAVE. Our choice of the Crank-Nicholson 
method is based on its recognized stability. The norm conservation is commonly used as a criterion 
to evaluate the reliability of the numerical calculations of the time evolution of wave functions. 



In References 



231 ] and 24]], this criterion is used to show analytically that the Crank-Nicholson 



method is unconditionally stable. Here, in order to evaluate the reliability of our algorithm, we 
have numerically calculated the norm of the wave packet for different times. The results thus 
obtained show that the norm is preserved. 

In fact, numerically one can only treat the tunneling from something process, where one gives a 
initial wave function with a given mean energy, very concentrated in a region next to x = 0. That 
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initial condition fixes an energy for the dust and the initial region where x may take values. Our 
choice for the initial wave function is the following normalized gaussian, 

^ (x , ) = ^^^y /4 xe(- 4 ^ 2 ), (46) 

where E is the dust energy. The wave- function \&(x,0) is normalized by demanding that the 
integral of |\&(3;,0)| 2 from to oo be equal to one and its mean energy be E. After one gives the 
initial wave function, one leaves it evolve following the appropriate Schrbdinger equation until it 
reaches infinity in the x direction. Numerically, one has to fix the infinity at a finite value, let us 
call it x maxc [. The general behavior of the solutions, when E is smaller than the maximum value 
of the potential barrier, is an everywhere well-defined, finite, normalized wave packet. Even in 
the limit of a vanishing scale factor. For small values of x the wave packet has great amplitudes 
and oscillates rapidly due to the interaction between the incident and reflected components. The 
transmitted component is an oscillatory wave packet that moves to the right and has a decreasing 
amplitude which goes to zero in the limit when x goes to infinity. As an example, we solve eq. 
(|44|) with A = 0.01. For this choice of A the potential barrier has its maximum value equal to 20. 
In order to see the tunneling process, we choose E = 19 for the initial wave function eq. (|46p . For 
that energy, we compute the points where it meets the potential barrier, the left (xu p ) and right 
(x r tp) turning points. They are, xi tp = 15.4085707 and x r t p = 26.933766. In the present case we 
fix x maxc i = 83.23582897, as the infinity in the x direction. In figure [IJ we show \^f(x, t max d)\ 2 for 
the values of A and E, given above, at the moment t max d = 21.5 when ^ reaches infinity. For 
more data on this particular case see Table U in the appendix. It is important to mention that the 
particular choice of numerical values for A and E, above and in the following examples, were made 
simply for a better visualization of the different properties of the system. 

Now, we compute the tunneling probability as the probability to find the scale factor of the 
universe to the right of the potential barrier. In the present situation, this definition is given by 
the following expression [is! ] . 

JZ \*(x,t max )\ 2 dx 
innt ^\^t max) \Hx' 140 
where, as we have mentioned above, the oo limit must be replaced by some suitably large value of 
x, for the sake of numerical calculations. 

Since, by normalization, the denominator of Eq. (|47j) is equal to the identity, TPi n t is effectively 
given by the numerator of Eq. (|47p . We consider, here, the dependence of TP on the energy E. 
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FIG. 1: \*ff(x, t max d)\ 2 = p, for A = 0.01, E = 19 at the moment t max d = 21.5 when ^ reaches infinity, located at 
x = 83.23582897. 



Therefore, we compute TPi n t for many different values of E for a fixed A. For all cases, we 
consider the situation where E is smaller than the maximum value of the potential barrier. From 
that numerical study we conclude that the tunneling probability grows with E for a fixed A. As 
an example, we consider 20 values of the dust energy for a fixed A = 0.01. For this choice of A the 
potential barrier has its maximum value equal to 20. In order to study the tunneling process, we 
fix the mean energies of the various ^(x,0)'s in Eq. (|46[) to be smaller than that value. Table U in 
the appendix shows the different energy values E, TPi nt , xu p and x r t p for each energy. 

Since TP grows with E it is more likely for the universe, described by the present model, to 
nucleate with the highest possible dust energy. Therefore, it is more probable that the classical 
evolution should start with the greatest possible value for the dust energy. 

Now, we may compare the above TP's for dust with the TP's for a radiative perfect fluid 
(a = 1/3). The TP's for a closed FRW model coupled to a radiative perfect fluid and a positive 



cosmological constant, described in Schutz's formalism, were first computed in Ref. [la ]. In order 
to compare the TP's between the two models we must fix values for A and E. Then, we must 
fix values for the infinity (x max ) in both models so that they be comparable. Here, we consider 
the distance from the point xq t where the potential vanishes (xor ^ 0) to x maxr , for the radiative 
model, as the reference. We compute the ratio x ma xr/xQ r . Let us call it A r . Then, we fix the 
Xmaxd, for the dust model, by multiplying xod, for this model, by A r . Finally, for fixed values of 
A, E and x max , we compute the TP's for both radiative and dust models, with the aid of eq. (I47p . 
We repeat these calculations for several different energy values E and conclude that the TP for 



14 



the dust is greater than the one for the radiation, for the same energy. As an example, we consider 
the same 20 values of the dust energy, given in table HI and repeat the calculations for the radiative 
model with fixed A = 0.01. For this choice of A the potential barrier, in the radiative model, has 
its maximum value equal to 225 and XQ r = 17.32050808. We choose x maxr = 30, in the radiative 
model. Therefore, A r = (30/17.32050808) = 1.732050807 and since x od = 48.05622828, for the 
dust model, x max d = 83.23582897. In table HU in the appendix, we can see the different energy 
values E, TPi nt , xi tp and x r t p for each energy, in the radiative model. Figure [2] shows the tunneling 
probability as a function of E in logarithmic scale, for both models. 
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FIG. 2: logTPint for different energies (E) for a fixed A = 0.01, for both radiative and dust models. 

IV. CONCLUSIONS. 

In the present work, we studied the quantum cosmology description of closed Friedmann- 
Robertson- Walker (FRW) models in the presence of a positive cosmological constant (A) and a 
generic perfect fluid. We worked in the Schutz's variational formalism. If one uses the scale 
factor and its canonically conjugated momentum as the phase space variables that describe the 
geometrical sector of these models, one obtains Wheeler-DeWitt equations with operator ordering 
ambiguities. In order to avoid such ambiguities and simplify the quantum treatment of the models, 
we introduced new phase space variables. We explicitly demonstrated that the transformation 
leading from the old set of variables to the new one is canonical. In the demonstration, carried out 
in the symplectic framework, we showed that the Dirac and Poisson brackets of the new canon- 
ical variables are the same. In order to show that the above canonical transformations simplify 



15 



the quantum treatment of those models, we considered a particular model. We applied it to the 
canonical quantization of a closed FRW model in the presence of a positive cosmological constant 
and dust. We used the Schutz's variational formalism and the appropriate canonical transforma- 
tion to re-write the Wheller-DeWitt equation as a Schrodinger equation, in the new variables, free 
from operator ordering ambiguities. In the present case, the Schrodinger equation had a potential 
barrier. We solved it numerically using the Crank-Nicholson scheme and determined the time 
evolution of the initial wave function. Then, using the wave function we computed the tunneling 
probability (TP) for the birth of an asymptotically DeSitter, inflationary universe, as a function of 
the mean energy E of the initial wave function. We observed that the TP grows with E. Finally, 
we compared the TP of the present model with the TP of another model already studied in the 
literature 18|]. There, the only difference is the presence of a radiative perfect fluid {a = 1/3), 
instead of dust. We observed that the TP for the dust is greater than the one for the radiation, 
for the same energy. 
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E 


TPint 


%it P 




19.000000 


0.0162028 


15.408570 


26.933766 


18.000000 


3.15309c-10 


13.130081 


29.397063 


17.000000 


3.98196c-19 


11.418328 


31.300238 


16.000000 


6.5559c-28 


10.004429 


32.912664 


15.000000 


2.27017e-36 


8.784433 


34.338844 


14.000000 


1.78018e-44 


7.705290 


35.632446 


13.000000 


3.07762e-52 


6.735782 


36.825393 


12.000000 


1.11065e-59 


5.855886 


37.938510 


11.000000 


7.86966c-67 


5.052047 


38.986278 


10.000000 


1.02929c- 73 


4.314784 


39.979254 


9.000000 


2.3364c-80 


3.637385 


40.925413 


8.000000 


8.63572e-87 


3.015151 


41.830959 


7.000000 


4.84268e-93 


2.444963 


42.700829 


6.000000 


3.77963e-99 


1.925061 


43.539032 


5.000000 


3.65203e-105 


1.454985 


44.348872 


4.000000 


3.65156e-lll 


1.035704 


45.133115 


3.000000 


2.47116e-117 


0.670037 


45.894101 


2.000000 


9.73233e-125 


0.363698 


46.633833 


1.000000 


1.2813e-130 


0.128371 


47.354041 


0.500000 


1.45689c-133 


0.045367 


47.707300 



TABLE I: The computed values of TPi n t, xu p and x rtp for 20 different values of E when A = 0.01, x max d = 
83.23582897 and t maxd = 21.5, for the dust model. 
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E 


TPint 


%it P 




19.000000 


6.4146e-296 


2.544209 


17.132630 


18.000000 


4.2118e-296 


2.474885 


17.142781 


17.000000 


2.68803e-296 


2.403736 


17.152902 


16.000000 


1.65734e-296 


2.330596 


17.162993 


15.000000 


9.79524e-297 


2.255268 


17.173054 


14.000000 


5.49363e-297 


2.177524 


17.183084 


13.000000 


2.8852e-297 


2.097094 


17.193086 


12.000000 


1. 39383c-297 


2.013655 


17.203058 


11.000000 


6.04368e-298 


1.926814 


17.213001 


10.000000 


2.27188e-298 


1.836087 


17.222914 


9.000000 


7.03829c-299 


1.740866 


17.232800 


8.000000 


1.664e-299 


1.640366 


17.242656 


7.000000 


2.65575e-300 


1.533548 


17.252485 


6.000000 


2.32358e-301 


1.418983 


17.262285 


5.000000 


7.59041c-303 


1.294616 


17.272057 


4.000000 


4.17488e-305 


1.157287 


17.281802 


3.000000 


5.36776e-309 


1.001676 


17.291519 


2.000000 


1.90162e-317 


0.817407 


17.301209 


1.000000 


5.92879c-322 


0.577672 


17.310872 


0.500000 


4.55125e-314 


0.408362 


17.315693 



TABLE II: The computed values of TPi nt , xit P and x rtp for 20 different values of E when A = 0.01, x maxr — 30.0 
and t m axr = 15.0, for the radiative model. 



